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O ■ Abstract 

(N 

We consider the minimum spanning tree (MST) problem under the restriction that for every vertex 
V, the edges of the tree that are adjacent to v satisfy a given family of constraints. A famous example 
thereof is the classical degree-bounded MST problem, where for every vertex v, a simple upper bound 
^^ \ on the degree is imposed. Iterative rounding/relaxation algorithms became the tool of choice for degree- 

constrained network design problems. A cornerstone for this development was the work of Singh and 
Lau ifTSl . who showed that for the degree-bounded MST problem, one can find a spanning tree violating 
C^ ' each degree bound by at most one unit and with cost at most the cost of an optimal solution that respects 

r^ \ the degree bounds. 

However, current iterative rounding approaches face several limits when dealing with more general 

Q \ degree constraints. In particular, when several constraints are imposed on the edges adjacent to a vertex 

V, as for example when a partition of the edges adjacent to v is given and only a fixed number of elements 

can be chosen out of each set of the partition, current approaches might violate each of the constraints 

^ ■ by a constant, instead of violating the whole family of constraints by at most a constant number of edges. 

On \ Furthermore, it is also not clear how previous iterative rounding approaches can be used for degree 

Cn ■ constraints where some edges are in a super-constant number of constraints. 

We extend iterative rounding/relaxation approaches both on a conceptual level as well as aspects 
involving their analysis to address these limitations. Based on these extensions, we present an algorithm 
r^ \ for the degree-constrained MST problem where for every vertex i\ the edges adjacent to v have to be 

independent in a given matroid. The algorithm returns a spanning tree of cost at most OPT such that 
for every vertex v, it suffices to remove at most 8 edges from the spanning tree to satisfy the matroidal 
degree constraint at v. 
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H ■ 1 Introduction 

a , 

Recently, much effort has been put on designing approximation algorithms for degree-constrained network 
design problems. This development was motivated by various applications as for example VLSI design, 
vehicle routing, and applications in communication networks |17][3l[l71- One of the most prominent and 
elementary problems here, which attracted lots of attention in recent years, are degree-constrained (MST) 
problems. 

In the most classical setting, known as the degree-bounded MST problem, the problem is to find a 
spanning tree T Q E of minimum cost in a graph G = {V,E) under the restriction that the degree of 
each vertex v with respect to T is at most some given value B^. Since checking feasibility of a degree- 
bounded MST problem is already NP-hard, interest arose in finding low-cost spanning trees that violate the 
given degree constraints slightly. A long chain of papers (see 171 [I2l [13] |4l HI and references therein) led 
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to algorithms with various trade-offs between cost of the spanning tree and violation of the degree bounds. 
In recent years, important progress was achieved for the degree-bounded MST problem, which also led to a 
variety of new techniques. Goemans t8j showed how to find a spanning tree violating each degree constraint 
by at most two units, and whose cost is bounded by the cost OPT of an optimum spanning tree that satisfies 
the degree constraints. Enhancing the iterative rounding framework introduced by Jain ^ with a relaxation 
step, Singh and Lau ifTSl obtained a stronger version of the above result, which is essentially best possible, 
where degree constraints are only violated by at most one unit. They work with an LP relaxation of the 
problem, and iteratively drop degree constraints from the LP that cannot be violated by more than one unit 
in later iterations. The adapted LP is then solved again to obtain a possibly sparser basic solution that allows 
for further degree relaxations. Edges not used in the current optimal solution to the LP are removed from the 
graph, and edges that have a weight of one are fixed, while updating degree bounds accordingly. A degree 
bound at a vertex v is removed whenever it is at most one unit lower than the current number of edges 
adjacent to v. 

We are interested in obtaining results of similar strength for more general degree bounds. Consider for 
example the following type of degree constraints: for every vertex v, a partition E^, . . . , E"^^ of the set 5{v) 
of edges adjacent to v is given, and within each set E'f of the partition, only a given number of edges can be 
chosen. The algorithm of Singh and Lau lITSl as well as the one of Goemans HI can easily be adapted to this 
setting. (In particular, the algorithm of Singh and Lau was even presented in this precise setting.) However, 
with both of these approaches, the constraint imposed by each set E"^ can be violated by a constant. We are 
interested in having at most a constant violation over all degree constraints at v, i.e., for every vertex v ^V , 
at most a constant number of edges have to be removed from the spanning tree to satisfy all constraints 
at V. Another more general example that will be useful to illustrate limits of cuiTcnt methods is obtained 
when imposing constraints for each vertex w on a laminar family on the edges adjacent to v, instead of only 
considering a partition. 

Adapting Goemans' algorithm to these stricter bounds on the degree violation seems to be difficult, since 
a crucial step of this algorithm is to cover the support E* of a basic solution to the natural LP relaxation 
by a constant number of spanning trees (for the degree-bounded MST problem, Goemans showed tBj that 
two spanning trees suffice). This result allows for orienting the edges in E* such that every vertex has at 
most a constant number of incoming arcs, at most one in each spanning tree. Dropping for every vertex 
all incoming arcs from its degree constraint then leads to a matroid intersection problem, whose solution 
violates each degree constraint by at most a constant. To be able to decompose E* into a constant number 
of spanning trees, one needs to show that for any subset of the vertices S* C F, only a linear number (in \S\) 
of edges have both endpoints in S. In the classical degree-bounded MST problem, this sparseness property 
follows from the fact that when considering only edges with both endpoints in S, there are at most a linear 
number (in |5|) of linearly independent and tight spanning tree constraints due to combinatorial uncrossing, 
and only a linear number of degree constraints within S. However, in more general settings as highlighted 
above, the number of degree constraints within S can be super-linear. 

Iterative relaxation looks more promising for a possible extension to generalized degree bounds. How- 
ever, current iterative rounding approaches face several limits when trying to adapt them. In particular, when 
dealing with the partition bounds as explained above, a simple adaptation of the relaxation rule, where for 
a vertex v all constraints at v would be dropped as soon as it is safe to do so due to a small support E* , 
risks to get stuck because there might be no vertex whose degree constraint can be relaxed. Furthermore, 
previous approaches (as used in ifTSl 13) to show that the support is spai^se fail in our setting because of a 
possible super-lineaiity of the total number of degree constraints. Additionally, previous iterative relaxation 
approaches crucially rely on the property that any edge is in at most a constant number of degree constraints 
to obtain violations that are bounded by a constant. However, this does not hold when dealing for example 
with degree constraints given by a laminar family. 

In this paper we show how to extend iterative relaxation approaches, both from a conceptual point of 



view as well as aspects involving their analysis, to tackle a wide class of MST problems with generalized 
degree bounds, namely when the degree bounds for every vertex are given by a matroid. In particular, this 
includes the partition bounds and the more general laminar bounds mentioned above. 

Our results. We present an iterative rounding/relaxation algorithm for finding a matroidal degree-bounded 
MST. The degree bounds are given as follows: for every vertex v, a matroid M^ = {6{v),Iy) over the 
ground set 5{v) is given with independent sets denoted by Xy C 2^^""^ The problem (without relaxed degree 
constraints) is to find a spanning tree T in G satisfying T n 5{v) G ly \/v € V, and minimizing a linear 
cost function c : E ^- IR+. We say that a given spanning tree T violates a degree constraint My by at most 
A; € N units, if it suffices to remove at most k edges R C d{v) fl T from T to satisfy the constraint My, 
i.e., {T\ R) n 5{v) € ly. Hence, the partition and laminar bounds mentioned above correspond to the case 
where all matroids My are partition or laminar matroids, respectively. We show the following. 

Theorem 1. There is an efficient algorithm for the matroidal degree-bounded MST problem that returns a 
spanning tree of cost at most the cost of an optimal solution, and violates each degree bound by at most 8 
units. 

To overcome problems faced by previous iterative relaxation approaches, we enhance the iterative relax- 
ation step, and exploit polyhedral structures to prove stronger sparseness results. The polytope used as a 
relaxation of the matroidal degree-bounded MST asks to find a point x € M^ in the spanning tree polytope 
such that for every vertex v G V, the restriction of x to 5{v) is in the matroid polytope P/vA, of My. 

To be able to always find possible relaxation steps, our iterative rounding procedure tries to achieve 
a somewhat weaker goal than previous approaches. The algorithm of Singh and Lau ifTSl relaxes degree 
constraints with the goal to approach the spanning tree polytope, which is integral. In our approach, the goal 
we pursue is to remove every edge {u, v} from at least one of the two degree constraints at u or v. As soon 
as no edge is part of both degree constraints at its endpoints, the problem is a matroid intersection problem, 
since all degree constraints together can be described by a single matroid over the support of the current 
LP solution. Thus, once we are in this situation, the current LP will be integral and no further rounding 
steps are needed. Hence, in our relaxation step, we try to find a vertex v such that we can remove all edges 
adjacent to v that are still in both degree constraints from the degree constraint at v. Edges adjacent to v that 
are only contained in the degree constraint at v will not be removed from the constraint My. Our approach 
has thus some similarities with Goemans' method, but instead of removing right at the start every edge from 
one degree constraint, we do this iteratively and hereby profit from additional sparseness that is obtained 
by solving the LP relaxation after each degree adaptation step. As we will see in Section |2l the way how 
we remove edges from a constraint is strictly speaking not a relaxation, and we therefore prefer to use the 
term degree adaptation instead of degree relaxation. The above degree adaptation step alone shows not to 
be sufficient for our approach, since one might still end up in a situation were no further degree adaptation 
can be performed because the graph is too dense. To obtain greater sparsity, we use a second type of degree 
adaptation, where for some vertex v we remove (almost) the full degree constraint at v if this cannot lead to 
a large violation of the degree constraint at v. 

The main step in the analysis is to prove that it is always possible to apply at least one of two suggested 
degree adaptations. A first step in this proof is to show that the support of a basic solution to the LP relaxation 
is sparse. We obtain sparsity by showing that if there are fc G N linearly independent and tight constraints 
(with respect to the current LP solution x) of the polytope Pm^, then x{5{v)) > k. Since summing x{6{v)) 
over all vertices is equal to 2{\V\ — 1), because x{E) = \V\ — I as x is in the spanning tree polytope, there 
are at most 2{\V\ — 1) linearly independent and tight degree constraints. 

The crucial part in the analysis is to show that vertices to which no further degree adaptation can be 
performed do not have very low degrees in average, implying that some of the other vertices are likely to 



have low degrees and therefore admit a degree constraint adaptation. To prove this property, we exploit the 
interplay between degree bounds and spanning tree constraints to show that any degree two node can either 
be treated separately and allows for reducing the problem, or implies a reduction in the maximum number 
of linearly independent and tight spanning tree constraints. 

Related work. The study of spanning trees with degree constraints can be traced back to Fiirer and 
Raghavachari |l7|, who presented an approximation algorithm for the degree-bounded Steiner Tree prob- 
lem which violates each degree bound by at most one, but does not consider costs. This result generated 
much interest in the study of degree-bounded network design problems, leading to numerous results and 
new techniques in recent years for a variety of problems, including degree-bounded arborescence problems, 
degree-bounded A:-edge-connected subgraphs, degree-bounded submodular flows, degree-bounded bases in 
matroids (see |[l6l[l7l[IIl[ni[T5l[l0ll2|6l[Dand references therein). 

Spanning tree problems with a somewhat different notion of generalized degree bounds have been con- 
sidered in Q and lUl. In these papers, the term "generalized degree bounds" is used as follows: given is a 
family of sets Ei, . . . , Ei^ C E, and the number of edges that can be chosen out of each set Ei is bounded 
by some given value Bi G N. In IH, using an iterative relaxation algorithm, whose analysis is based on 
a fractional token counting argument, the authors show how to efficiently obtain a spanning tree of cost at 
most OPT and violating each degree bound by at most maxgg^; |{i € [/c] | e E £^j}|, the maximum coverage 
of any edge by the sets E^. In fill, a new iterative rounding approach was presented for the problem when the 
sets El, . . . ,Ek correspond to the edges E-i = 5{Ci) of a family of cuts Ci C V for i G [k] that is laminar. 
Contrary to previous settings where iterative rounding approaches were applied, here, it is possible that an 
edge lies in a super-constant number of degree constraints. At each iteration, the algorithm reduces the num- 
ber of degree constraints by a constant factor, replacing some constraints with new ones if necessary. This 
is done in such a way that degree constraints are violated by at most a constant in every iteration, leading to 
a spanning tree of cost at most OPT, that violates each degree constraint by at most 0(log(|F|)). 

Organization. In Section [2] we present our algorithm for the matroidal degree-bounded MST problem. 
The analysis of the algorithm is presented in Section [3] 

2 The algorithm 

Since during the execution of our algorithm the underlying graph will be modified, we denote hy H = 
{W, F) the current state of the graph, whereas G = {V, E) always denotes the original graph. For brevity, 
terminology and notation is with respect to the current graph H when not specified further. To distinguish 
between initial degree constraints and current degree constraints, we denote by N^ the current constraints 
for w G W — which will as well be of matroidal type — whereas M^ denotes the initial degree constraints 
at w € y. The vertices of H are called nodes since they might contain several vertices of G due to edge 
contractions. 

The algorithm starts with H = G and A^^, = M^ for v ^ V, and the LP relaxation we use is the 
following, 
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where Pgt denotes the spanning tree polytope of H, P^^^ denotes the matroid polytope that corresponds 

to N.,,, and xL, , 
correspond to 6{w 



to A^u,, and xL, . denotes the vector obtained from x G M^ by considering only the components that 



Algorithm for Matroidal Degree-Bounded Minimum Spanning Trees 

1. Initialization: H = {W, F) ^ G = {V, E), N^ ^ My for v eV. 

2. While |H^| > Ido 

a) Determine basic optimal solution x to (LPl). Delete all edges f ^ F with x{f) = 0. 

b) Contract all edges f e F with x{f) = 1. 

c) Fix a maximal family of linearly independent and tight spanning tree constraints. 

d) Type A degree adaptation: for each node w (^ W such that the set of all edges U C 5{w) 
that are still in both degree constraints is non-empty and satisfies \U\ — x{U) < 4, 
remove U from the degree constraint Ny,. 

e) Type B degree adaptation: for each node w ^W such that the set of all edges U C 5{w) 
contained in the degree constraint A^^ but not adjacent to a node in Q is non-empty and 
satisfies \U\ — x{U) < 4, remove U from the degree constraint Nw 

3. Return all contracted edges. 



There is a set of nodes Q = Q{H,x) C W that has a special role in our algorithm due to its relation 
with tight spanning tree constraints. The node set Q is defined and used after having contracted edges of 
weight one. Hence, assume that H does not contain any edge f & F with x{f) = 1. Then Q is defined as 
follows: we start with Q = and as long as there is a node w € W such that x{6{w) n F[W \ Q]) = 1, 
where F[W \ Q] is the set of all edges with both endpoints in W\Q, we add w to Q. One can easily observe 
that Q does not dependent on the order in which nodes aie added to io As we will see later, edges adjacent 
to these nodes can often be ignored from degree constraints due to strong restrictions that are imposed by 
the spanning tree constraints. 

The box on top of the page gives a description of our algorithm, omitting details of how to deal with 
the matroidal degree bounds when removing or contracting edges. We discuss these missing points in the 
following. 

Notice, that a basic solution to {LPl) can be determined in polynomial time by the ellipsoid method, 
even if the involved matroids are only accessible trough an independence oracle. Depending on the matroidal 
degree bounds involved, (LPl) can be solved more efficiently by using a polynomially-sized extended 
formulation. 

A tight spanning tree constraint, as considered in step (l2c] ). corresponds to a set L C W^, L ^ such 
that x{F[L]) = \L\ — 1. Fixing a tight spanning tree constraint means that this constraint has to be fulfilled 
with equality in all linear programs of type (LPl) solved in future iterations. It is well-known that if 
supp{x) = F, then any maximal family of linearly independent and tight spanning tree constraints with 
respect to x defines the minimal face of the spanning tree polytope on which x lies (see e.g. H). Hence, 
due to step (l2cl ). we have that if the LP solution at some iteration of the algorithm is on a given face of the 
spanning tree polytope, then all future solutions to {LPl) will be as well on this face. 

Fixing tight spanning tree constraints shows to be useful since they often imply strong conditions on 
the edges, which can be exploited when having to make sure that degree constraints ai^e not violated too 
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much. In particulai^, consider a node w ^ Q which, in the iterative construction of Q, could have been 
added as the first node, i.e., x{5{w)) = 1. When fixing tight spanning tree constraints, one can observe that 
any spanning tree satisfying those tight constraints with equality contains precisely one edge adjacent to w. 
Furthermore, the fixing of tight spanning tree constraints guai^antees that a node w G Q will stay in Q in 
later iterations until an edge adjacent to w is contracted. Hence, all edges being in some iteration adjacent to 
a node w £ Q, will be adjacent to a node in Q in all later iterations until they are either deleted or contracted. 
This property is important in our approach since a type B degree adaptation ignores edges adjacent to Q, 
and we want to make sure that an edge which is once ignored will never be considered during a later type B 
degree adaptation. 

Contracting and removing edges. To fill in the remaining details of our algorithm, it remains to discuss 
how edges are contracted and removed. Throughout the algorithm, any degree constraint N^ of a node w 
containing the vertices vi, . . . ,Vk € V can always be written as a disjoint union of matroidal constraints 
Ny-^ , • • • , Ny^ , where Ny. corresponds to the "remaining" degree bound at Vi and is a matroid over the edges 
6{w) that are adjacent to Vi. Whenever an edge / = {wi,W2} of weight one is contracted in step (l2bl ) of 
the algorithm to obtain a new node w, the new degree constraint A^^ at w is obtained by taking a disjoint 
union of the matroids Ny^^/ f and N^^/ f, where Nyj^/ f and N^^j j correspond to the matroids obtained 
from Nyj^ and Ny,^, respectively, by contracting /. This operation simply translates the degree constraints 
on w\ and W2 to the merged node w. The property that a degree bound on w is a disjoint union of degree 
bounds of the vertices represented by w, is clearly maintained by this contraction. 

As highlighted in the box, we adapt constraints by removing for some node w € VF a set of edges 
U C b{w) from the constraint Nw When removing U from Ny,, we construct a new degree constraint given 
by a matroid N'y, over the elements b{w) such that the following properties hold. 

Property 2. 

i) N^ is a disjoint union of matroidal constraints Ny_^ , • • • , N' corresponding to vertices contained in w, 

ii) edges ofU are free elements of N[y, i.e., if I is independent in N^ then I L)U is independent in N'^,, 

Hi) any independent set ofNy. can be transformed into one ofNy. by removing at most \\U\ — x(U)~\ edges, 

iv) the previous LP solution x is still feasible with respect to N'y,, i.e., x\r, . € Pj 
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Any removal operation satisfying the above properties can be used in our algorithm. Before presenting 
such a removal operation, we first mention a few important points. To avoid confusion, we want to highlight 
that removing U from N^ does not simply correspond to deleting the elements U from the matroid N^. For 
any edge / € 5{w) that is free in Ny^, we say that / is not contained in the degree constraint Nyj, and it is 
contained otherwise. When all edges adjacent to a given node w are not contained in its degree constraint, 
which corresponds to N^ being a free matroid, we say that the node w has no degree constraint. 

We now discuss how to remove a set of edges U C 8{w) from Ny, to obtain an adapted degree bound 
N'y, satisfying Property |2] Let vi, . . . ,Vp € F be all vertices contained in the node w, and we consider the 
decomposition of Ny, into a disjoint union of matroids Ny^ , • • • , Ny^, , where Ny^ for i € [k] corresponds 
to the "remaining" degree bound at Vi. To remove U from Ny,, we adapt each matroid A'^;,- as follows to 
obtain a new matroid A'^^ . Let Si be the ground set of Ny., i.e., all edges in 5{w) being adjacent to Vi. Let 
-^1 = ('S'i,2^i) be the matroid with independent sets 

ii = {iQSir\U\\i\< \Si r\u\- \x{Si n u)\}. 

Hence, Mi is a special case of a partition matroid. Let M2 = Mi V A'^- be the union of the matroids Mi 
and Ny-, and let M3 = M2/{Si n U) be the matroid obtained from M2 by contracting Si fl U. The degree 
bound A'^^, is obtained by a disjoint union of M3 and a free matroid over the elements in Si n U. The new 
degree constraint A^^,, that results by removing U from N^;, is given by the disjoint union of the matroids 

A^' A^' 



Lemma 3. TJie above procedure to remove elements from a degree constraint satisfies Property^ 

Proof. By construction, when removing a set [/ C 5{w) from a degree bound N^, which can be written 
as a disjoint unions of N^^ , • • • , N^^ , a matroidal bound A'^^ is determined which is a disjoint union of 
N^^ , . . . , AT^^ . Hence point (0) of Property [2] holds. 

Let Si be the ground set of the matroids N'j;.,Ny^ for i G [k]. Since A'^^ is a disjoint union of 
N^^, . . . , Ny^, it suffices for point ^ to prove that if /' is independent in N^. then P U {Si n U) is in- 
dependent in NIj_. This follows since N^_ was obtained by a disjoint union of the matroid M3, as defined 
above, and a free matroid over Si n U. 

For point (Iml ). consider an independent set /' in N^.. Since all edges in U H Si are free in N'^., we 
can assume {U H Si) C /'. Consider how the matroid N^_ was constructed by the use of the matroids 
Ml, M2, M3. We start by observing that U H Si is an independent set in M2 = Mi V A^^^,-. Let r^ be the rank 
function of A^^,-, and r2 be the rank function of M2. Since x G -P/v„. , we have that ri{Si nU) > x(S'j PI C/). 
Furthermore, since M2 = Mi V A'^, and any jS'j n L'^l — [x(5i D U)\ elements of Si nU are independent 
in Ml, we have 

r2iSi nU) = min{\Si D U\,ri{Si nu) + \Si nu\- [xiSi nu)\} = \Si n u\, 

showing independence of 5/ CiU in M2. Because A'^^ was obtained by a disjoint union of the matroid M3 
and a free matroid over the elements Si H U, we can write P = I^U {Si fl U) with I3 independent in M3. 
Furthermore, as M3 = M2/{Si fl U) and Si H U is independent in M2, the set /' is independent in M2. 
As M2 = Ml V Ny^, we have /' = /i U /, with /i independent in Mi and / independent in A'^,^. Since 
Ml is a matroid of rank \Si n f/| — [x{Si n U)\, we have that / is obtained from /' by removing at most 
\h\ < \Si nU\- [x{Si n U)\ < \U\ - [x{U)\ elements as desired. 

Let Xj = xL for i e [A;]. To show point (jlvll, it suffices to prove that Xj € P/y/ \/i € [A;], since N' is a 
disjoint union of N!^_^ , • • • , A^^^^.. Let Zi G [0, 1]"^' be given by 

[xtU) i!f€S.\U^ 

Observe that Zj — Xj G Pmj because the support of Zj — Xj is a subset of Si nU, \\zi — Xi\\i = \SifMJ\ — 
x{Si n U) and any \Si nU\ — [x{Si n U)\ elements of Si nU are independent in Mi. Hence Zi G Pm2. 
since M2 = Mi V A'^y., Xj G P/v„, and Zj — Xj G Pmi- As M3 = M2/{Si fl C/), we have that the restriction 
of Zi on Si \ U, which is equal to Xj| axri, is in Pm^- Since A'^^',. is the union of M-^ and a free matroid over 
Si n U, this finally implies that Xj G P/v . 

D 

3 Analysis of the algorithm 

Lemma 4. During the execution of the algorithm, for every vertex v ^V, at most one constraint adaptation 
of type A and one of type B is performed that removes edges of 5 {v) from degree constraints containing v. 

Proof. When a type A degree adaptation is applied to a node w ^ W that contains v, no further type A 
degree adaptation can remove any edges in 6{v) n F from the constraint containing v, since those edges are 
not anymore contained in both degree constraints at their endpoints. 

Similarly, when a type B degree adaptation is applied to a node w that contains v, all edges in 6{w) n 
F[W \ Q] are removed from the degree constraint at w and thus cannot be removed again at a later type 
B degree adaptation. Hence, the only possibility to remove further edges adjacent to f in a later type B 



degree adaptation is that some edge f ^ F which was — at some iteration of the algorithm — not considered 
for a possible removal by a type B adaptation because of being adjacent to a node in Q, can be removed 
by a type B adaptation at a later stage. However as already discussed, since we fix all tight spanning tree 
constraints, an edge that is adjacent to a node w £ Qin some iteration, will remain so until it is either deleted 
or contracted in step (l2al ) or (l2bl ) of the algorithm. Hence, this "bad constellation" can never occur. D 

We exploit that our removal operation satisfies point (Iml ) of Property |2] to bound the maximum possible 
degree violation. In particular, for each vertex v G V, every time edges U with U n 5{v) 7^ are removed 
from the current degree constraint N,,^ at the node w that contains v, the degree constraint at v can be 
violated at most by an additional \\U\ — x{U)'] units. Since we only perform degree adaptations for sets U 
with \U\ — x{U) < 4, and Lemma |4] guarantees that at most two adaptations are performed that involve the 
degree constraint at t", we obtain the following result. 

Corollary 5. If the algorithm terminates, then the returned tree violates each degree constraint by at most 
8 units. 

A main step for proving that we can always apply one of the two suggested degree adaptations, is to 
prove that a basic solution to {LPl) is sufficiently sparse. A first important building block for proving 
sparsity is the following result. 

Lemma 6. Let x be any solution to {LPl) whose support equals F. Then for every node w € W, the 
maximum number of linearly independent constraints of the matroid polytope Pn^, that are tight with respect 
to X, is bounded by x{6{w)). 

Proof. Let C C 2^^'"^ be a family with a maximum number of sets that correspond to linearly independent 
constraints of the matroid polytope i-jy^^, that are tight with respect to x. By standard uncrossing arguments, 
C can be chosen to be a chain, i.e. C = {Ci, . . . , Cp} with Ci C C2 C • • • C Cp (see |3 for more details). 
We have to show that p < x{6{w)). Let r be the rank function of N^,. Define Co = and for i € [p] let 
Ri = Ci\ Ci-i. Since C is a family of tight constraints, we have 

x{Ili) = r{Ci) - r{a.i) yie[p]. (1) 

Because Ri C supp{x), the left-hand side of ([T]) is strictly lai^ger than zero. Furthermore, the right-hand 
side is integral and must therefore be at least one. Hence x{Ri) > 1 f or i G \p\, which implies x{5{v)) > 

TIl=l<R^)>p■ n 

Notice, that the above lemma implies that a basic solution x to {LPl) has a support of size at most 
3(|T4^| — 1), because of the following. We can assume that all edges that are not in the support of x are 
deleted from the graph. Due to Lemma[6l at most Yliwe.w x{5{w)) linearly independent constraints of the 
polytopes {Pnu, I ^ € W} can be tight with respect to x, and since x is in the spanning tree polytope of 
H, this bound equals X^^gvK^('^(^)) ~ 2(|VF| — 1). Furthermore, at most \W\ — 1 lineai^ly independent 
constraints of Pst ai^e tight with respect to x due to uncrossing. This shows in particular that in the first 
iteration of the algorithm, we can find a node w ^W io which a type A degree constraint adaptation can be 
applied, because 

Y^ {\6{w)\ - x{5{w)) = 2\F\ - 2{\W\ - 1) < 4(|VF| - 1), 

■w£W 

and hence there must be a node w ^W with \6{w) \ — x{5{w)) < 4. 

However, in later iterations, the above reasoning alone is not anymore sufficient because many vertices 
do not have degree constraints anymore. Still, by assuming that no type B constraint adaptation is possible. 



and using several ideas to obtain stronger sparsity, we show that the above approach of finding a good vertex 
for a type A degree adaptation by an averaging argument can be extended to a general iteration. 

For the rest of this section, we consider an iteration of the algorithm at step (l2dl) with a current basic 
solution X to (I/Pl), and assume that \W\ > 1, and that no type B degree adaptation can be applieco We 
then show that there is a type A constraint adaptation that can be performed under these assumptions. This 
implies that our algorithm never gets stuck, and hence proves its correctness. 

Since we often deal with the spare 1 — x(/) of an edge / S F, we use the notation z = 1 — x. 
Furthermore, we partition F into the sets F2,Fi and Fq of edges that are contained in 2, 1 and de- 
gree constraints, respectively. Hence, at the first iteration we have F2 = F. Our goal is to show that 
J2wew z{S{w) n F2) = 2z{F2) < 4|y I, where Y C W is the set of all nodes w with 5{w) n F2 / 0. By 
an averaging argument this then implies that there is at least one node w G Y to which a type A constraint 
adaptation can be applied. Notice that the set F2 cannot be empty (and hence also y 7^ 0): if F2 = 0, then 
the current LPl con^esponds to a matroid intersection problem since every edge is contained in at most one 
degree constraints, and hence all degree constraints form together a single matroid over F; in this case LPl 
is integral and a full spanning tree would have been contracted after step (l2bl ). which leads to \W\ = 1 and 
contradicts our assumption \W\ > 1. 

Lemma 7. Let Cbe a maximum family of linearly independent spanning tree constraints that are tight with 
respect to x. Then 

2z{F2) < 2\C\ + 2{\W\ - 1) - 2(|Fo| + |Fi|) - 2x(Fo). 

Proof. We can rewrite 22; (F2 ) as follows by using the fact that x (F) = \W\ — 1 (because x is in the spanning 
tree poly tope of H). 

2z{F2) = 2z{F) - 2z{Fo) - 2z{Fi) 

= 2{\F\-x{F))-2z{Fo)-2z{Fi) (2) 

= 2\F\ - 2{\W\ - 1) - 2z{Fo) - 2z{Fi) 

Using classical arguments we can bound the size of the support of x, which is by assumption equal to \F\, 
by the number of linearly independent tight constraints from the spanning tree polytope and the degree 
polytopes P/v^ for w S W . In particular x is uniquely defined by the tight spanning tree constraints C 
completed with some set V of linearly independent degree constraints, and we have |F| = j£j + \V\. The 
degree constraints V can be partitioned into V^ for w G W , where V^ are linearly independent constraints 
of the matroid polytope -P/v„. By Lemma|6l \'Du,\ is bounded by the sum of x over all edges in 5{w) that are 
contained in the degree constraint at w. When summing these bounds up over all w G W, each edge in F2 
is counted exactly twice, and each edge in Fi exactly once. Hence, 

\V\ < 2x{F2) + x{Fi) = 2x{F) - x{Fi) - 2x{Fq) = 2{\W\ - 1) - x{Fi) - 2x{Fq). 

Using |F| = \C\ + |Pj and the above bound, we obtain from ^ 

2z{F2) < 2\C\ + 2{\W\ - 1) - 2 (z(Fo) + ^(i^i) + 2x{Fo) + x{Fi)) 
= 2\C\ + 2(|Ty| - 1) - 2(|Fo| + li^il) - 2x{Fo), 

where the last inequality follows from z{U) + x{U) = \U\ for any U (^ F. D 

The size of a family C of linearly independent tight spanning tree constraints can easily be bounded by 
I VF| — 1 using the fact that one can assume C to be laminai^ by standard uncrossing arguments (and C contains 



^Notice that the assumption \W\ > 1 is not redundant. Whereas we know that at the beginning of the iteration | PF | > 1 did 
hold, this could have changed after contracting edges in step l l2bt . 



no singleton sets). However, this result shows not to be strong enough for our purposes. To strengthen this 
bound we exploit the fact that if L contains close to | VF| — 1 sets, then there are many nodes w ^W that are 
"sandwiched" between two sets of £, i.e., there are two sets Li, L2 € £ with L2 = Li U {w\, which in turn 
implies x(b{w) PI E\L'2\) = 1. Notice that for any degree two node w which is not in Q, we have x{U) 7^ 1 
for all U C 6{w). Hence, such a node cannot be "sandwiched" between two tight spanning tree constraints, 
and we expect that the more such nodes we have, the smaller is \C\. The following result quantifies this 
observation. It is stated in the general context of a spanning tree polytope of a general connected graph (not 
being linked to our degree-constrained problem). 

Lemma 8. Let y be a point in the spanning tree polytope for a given graph G = (V, E), and let 

S{G,y) = {v(^V\ \5{v)\ = 2, y{U) / 1 VC/ C 5{v)}. 

Let C Q2^ be any linearly independent family of spanning tree constraints that are tight with respect to y. 
Then 

\C\<\V\-l- \\S{G,y)\ 

Proof. To simplify notation let S = S{G, y). By standai^d uncrossing ai^guments (see for example [H), we 
can assume that C is laminar. We first consider the case that there is a set L € £ with L C S. Let L be a 
minimal set in C with this property. Since L is a tight spanning tree constraint, we have that y| Rrri is in the 
spanning tree polytope of G[L], and hence y{6{v) fl E[L]) > 1 for w E L. As L C 5, we have \5{v)\ = 2 
and y{e) < Ifor V £ L and e G S{v). This implies that every vertex in L must have both of its neighbors in 
L to satisfy y{6{v) H E[L]) > 1. Since G is connected, as we assumed that there is a point in the spanning 
tree polytope of G, we must have L = V = S. Furthermore \V\ > 3, because vertices in L have degree two. 
Hence the claim trivially follows since \C\ = 1. 

Now assume that there is no set L G £ with L C 5. We show that there exists a set i? C 5 of size at 
least \R\ > ^\S\, such that the laminai^ family Cr = {L\ R \ L £ C} over the elements V \ R satisfies the 
following: 
i) Cr has no singleton sets, 

ii) \Cji\ = \C\, i.e., any two sets Li, L2 £ C with Li C L2, satisfy L2 \ /^i 2 R. 
Notice that this will imply the claim since \Cr\ < \V \ R\ — I, because JZr is laminar without singleton 
sets, and hence j£j = \Cr\ < |F \ -R| — 1 < jV^| — 1 — ^\S\. It remains to define the set R with the desired 
properties. For L £ C, let Vl Q Lhe. all vertices in L that are not contained in any set P G £ with P C L. 
For each set L G £, include an arbitrary set of [^15 fl Vl\] elements of 5 n Vl in R. Since the sets Vl 
for L G £ are a partition of all vertices V, we clearly have \R\ > ^\S\. Furthermore R satisfies the desired 
properties as we show below. 

i) Assume by sake of contradiction that Cr contains a singleton set, i.e., there is a set L G C with 
\L\R\ = 1. We can assume that L is a minimal set in C By assumption we have L ^ S, and since R'^ S, 
the element in L \ i? is not in S. Hence, R contains all elements LnS, which is only possible if |L n 5| = 1 
and therefore \L\ = 2. However, this implies that there must be an edges of weight one between the two 
vertices in L, which contradicts the fact that one of those vertices is in S. 

ii) Assume by contradiction that there are two sets Li, L2 G £ with Li C L2 that satisfy -L2 \ ^1 ^ R- 
We can choose Li and L2 such that there is no set L G £ with Li C L C L2. By choice of R, this can only 
happen if L2 \ Li contains exactly one vertex v £ S. This implies y{S{v) n E[L2]) = 1, which contradicts 
the fact that v € S. D 

Lemma[8]can easily be generalized to the subgraph of a given graph G obtained by deleting the vertices 
Q{G, y). This form of the lemma is more useful for our analysis because of our special treatment of vertices 
inQ. 
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Lemma 9. Let y be a point in the spanning tree polytope of a given graph G = (V^ E) with y{e) 7^ 1 Ve € -E, 
let G' = G[V \ Q{G, y)], and let y' be the projection of y to the edges in G'. Let C be any linearly 
independent family of spanning tree constraints ofG that are tight with respect to y. Then 



\C\ <\V\-l 



\\S{G\y') 



Proof By standard uncrossing arguments, we can assume that £ is a maximal laminar family of tight 
spanning tree constraints. We prove the result by induction on the number of elements in Q = Q{G, y). If 
Q = 0, then the result follows from Lemma[8] Let g € Q be a possible first element added to Q during the 
iterative construction of Q, i.e., y{6{q)) = 1. This implies that V \ {q} is a tight spanning tree constraint. 
Let H = G[V \ {q}], yu = y\E[v\{q}] ^"'^ Qh = Q{H,yH)- Since Qh = Q\ {q}, we can apply the 
induction hypothesis to the graph H to obtain that any maximal family Ch of linearly independent tight 
spanning tree constraints in H with respect to yn satisfies \Ch\ ^ IV \ {q}\ — 1 — [||5(G',y')|J . The 
claim follows by observing that C = Ch U {V} is a maximal family of tight spanning tree constraints in G, 
and hence 

|£| = |£i:^| + l < |y|-l ^' 



i^\S{G',y') 



D 



Combining Lemma |9] with Lemma |7] we obtain the following bound, where we use S = S{H\W \ 
Q],x|^r^>„,) to simplify the notation. To get rid of the rounding on i|5| we use 2[i|S'|J > \S\ — 1. 

Corollary 10. 

2z{F2) < 4{\W\ - 1) - 2(|Fo| + iFil) - 2x(Fo) - \S\ + 1. 

The following lemma implies the coiTcctness of our algorithm. We recall that Y C W is the set of all 
nodes w eW such that 5{w) n F2 / 0. 

Lemma 11. There is a node w (zY such that a type A constraint adaptation can be applied to w. 

Proof l^tiY = W\Y. We will prove that 

4|F| < 2(|Foj + |Fi|) + 2x(Fo) + \S\. (3) 

Together with Corollary [TO]this then implies 2z{F2) < 4|y | — 3, which in turn implies by an averaging 
argument that there is at least one node in Y to which a type A constraint adaptation can be applied. To 
prove ^ we apply a fractional token counting ai^gument: we show that if we interpret the right-hand side 
of (l3]l as a (fractional) amount of tokens, then we can assign those tokens to the vertices in Y such that each 
vertex in Y gets at least 4 tokens. 

We think of the tokens corresponding to 2 ( | Fq | + 1 Fi | ) + 2a; (Fq ) as residing at the endpoints of the edges 
in Fq U Fi. Each edge / G Fq gets 2 + 2x(/) tokens, 1 + x{f) at each endpoint. Each edge / € Fi gets 
1 + x{f) tokens at the endpoint which does not contain / in its degree constraint, and 1 — x{f) tokens at the 
other endpoint. The tokens assigned to the endpoints of the edges thus sum up to 2(|Fo| + |Fi|) + 2x(Fo). 

We start by assigning tokens to vertices in Q. By definition of the vertices in Q, we can order the 
elements mQ = {qi, . . . , qp} such that for i € [p], we have x{FqJ = 1 where Fq^ = {{qi,v} S F j -y G 
W \ {qi, . . . , qi-i}}. Since x(FqJ = 1 and no edge f G F satisfies x{f) = 1 (such an edge would have 
been contracted), we have \Fq. \ > 2. Each vertex q-i £ Q gets all the tokens at both endpoints of the edges 
in Fq.. Since \Fq^ | > 2, gj receives indeed at least four tokens. 
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Let H' = {W' , F') = H\W\Q]he the induced subgraph over the vertices W \Q , and let x' = x\p,. 
Notice that x' is in the spanning tree polytope of H' since the set of edges U (^ F that have at least one 
endpoint in Q satisfy x{U) = \U\, and hence x'{F') = |F'| — 1. 

The remaining tokens are allocated as follows. Each node w ^Y r\ W gets for every edge / G 5^' (w), 
the tokens of / at the endpoint at w. Furthermore, every node in S gets an additional token from the term 
|5|. 

The attributed tokens clearly do not exceed the right-hand side of (l3]l. It remains to show that each node 
w £ Y n W gets at least 4 tokens. We distinguish the following three cases: (i) w (^ S, (ii) w ^ S and 
none of the edges 6h' (w) is contained in the degree constraint at v, and (iii) w ^ S and at least one edge of 
6h' (w) is contained in the degree constraint at w. Notice that the vertices considered in case (i) are precisely 
all vertices in H' of degree two, because if there was a degree two vertex w G W \ S, then w would have 
been included in Q. Hence, all vertices considered in case (ii) or case (iii) have degree at least 3 in H'. 

Case (i): w G S. Because \5h'{w)\ = 2, we have that both edges in 5h'{w) are not contained in the 
degree constraint at w, since otherwise a type B degree adaptation could have been performed at w. Hence, w 
receives 2+j;(/i)+2;(/2) tokens from those two edges plus one token from |5|, resulting in 3+a;(/i)+2;(/2) 
tokens. Since x' is in the spanning tree polytope of H', we have x(/i) + x(/2) = x{5h'{w)) > 1, and thus 
w receives at least 4 tokens. 

Case (ii): w ^ S and none of the edges 5h'{w) is contained in the degree constraint at w. The total 
number of tokens received by w thus equals \5h'{w)\ + x{6h'{w)) > 3 + x{5h'{w)), since \5hi{w)\ > 
3. The claim follows again by observing that x' is in the spanning tree polytope of H', which implies 

X^dH'iw)) > 1. 

Case (iii): w ^ S and at least one edge of 6h' (w) is contained in the degree constraint at w. Let U 
be the set of all edges in 6h' [w] that are contained in the degree constraint at w. Since no type B degree 
adaptation can be performed at w, we have \U\ — x{U) > 4. However, \U\ — x{U) is exactly the number of 
tokens that w receives from the edges in U. Hence, at least 4 tokens are assigned to w. D 
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